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A study of mixing and indirect CP violation in mesons through the determination 
of the parameters ycp and is presented. The parameter ycp is the deviation from 
unity of the ratio of effective lifetimes measured in decays to the CP eigenstate 
K~^K~ with respect to decays to the Cabibbo favoured mode K~'k~^ . The result 
measured using data collected by LHCb in 2010, corresponding to an integrated 
luminosity of 29pb~^, is 

ycP = (5.5 ± 6.3stat ± 4.1syst) x lO'^. 

The parameter is the asymmetry of effective lifetimes measured in decays of 
and mesons to K'^K" . The result is 

= (-5.9 ± 5.9stat ± 2.1syst) x lO-^. 
A data-driven technique is used to correct for lifetime-biasing effects. 
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1 Introduction 



Mixing of neutral mesons has only recently been established [H El [3] and first evidence 
for CP violation in the charm sector has just been seen by LHCb [1]. In this work the 
mixing and CP violation parameters ycp and A-p in the decays of neutral mesons into 
two charged hadrons are studied. Both quantities are measured here for the first time at 
a hadron collider. The observable ycp is the deviation from unity of the ratio of inverse 
effective lifetimes in the decay modes — )■ K^K~ and — K^ti^ 

_ t{D^ -> K+K~) 
Vcp = — 1) (1) 

where effective lifetime refers to the value measured using a single exponential model. All 
decays implicitly include their charge conjugate modes, unless explicitly stated otherwise. 
Similarly, is given by the asymmetry of inverse effective lifetimes as 

Ay = A= -. (2) 

r(L)0 ^ K+K-) + r(L)0 ^ K+K~) 

The neutral mass eigenstates |-Di,2) with masses mi^2 and widths ri^2 can be 
expressed as linear combinations of the flavour eigenstates as |-Di,2) = p\D'^) ±g|D°) with 
complex coefficients p and q satisfying |pp + \q\^ = 1. The average mass and width are 
defined as m = (mi + m2)/2 and F = (Fi + F2)/2; the mass and width difference are used 
to define the mixing parameters x = (m2 — mi)/F and y = (F2 — Fi)/(2F). The phase 
convention is chosen such that CP|Z)°) = —\D^) and CP\D^) = —\D'^) which leads, in 
the case of no CP violation [p = q), to \Di) being the CP odd and ID2) the CP even 
eigenstate, respectively. 



The parameter 



A 



A 



(3) 



contains the amplitude Af (Af) of (D^) decays to the CP eigenstate / with eigenvalue 
rjcp- The mixing parameters x and y are known to be at the level of 10~^ while both 
the phase and the deviation of the magnitude from unity of Aj are experimentally only 
constrained to about 0.2 [S]. The direct CP violation, i.e. the difference in the rates of 
and decays, is constrained to the level of 10~^ and has recently been measured by 
LHCb [1]. Introducing \q/p\^'^ ^ 1 ± Am and lAj/A/l^'^^ ^ 1 ± Ad, with the assumption 
that Am and Ad are small, and neglecting terms below 10~^ according to the experimental 
constraints, one obtains according to Ref. [B] 

ycp ~ + ^ Am^^ y cos - ^Amo; siu 0. (4) 

In the limit of no CP violation ycp is equal to y and hence becomes a pure mixing 
parameter. However, once precise measurements of y and ycp are made, any difference 
between y and ycp would be a sign of CP violation. 
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Previous measurements of ycp have been performed by BaBar and Belle. The results 
are ycp = (11.6 ± 2.2 ± 1.8) x lO'^ [7| for BaBar and ycp = (13.1 ± 3.2 ± 2.5) x lO'^ |2] 
for Belle. They are consistent with the world average oi y = (7.5 ±1.2) x 10~^ [5]. 

The study of the lifetime asymmetry of and mesons decaying into the singly 
Cabibbo-suppressed final state K^K~ can reveal indirect CP violation in the charm 
sector. The measurement can be expressed in terms of the quantity A^. Using the same 
expansion as for ycp leads to 



Ar 



i(A. 



+ Af{)y cos — X sin < 



ycp 



- (Am + Ad)y cos — X sin I 



(5) 



Despite this measurement being described in most literature as a determination of indirect 
CP violation it is apparent that direct CP violation at the level of 10~^ can have a 
contribution to at the level of 10~^. Therefore precise measurements of both time- 
dependent and time-integrated asymmetries are necessary to reveal the nature of CP 
violating effects in the system. 

The measurement of Ar requires tagging the flavour of the at production, which 
will be discussed in the following section. Previous measurements of A^ were performed 
by Belle and BaBar leading to Ar = (0.1 ± 3.0 ± 1.5) x lO'^ |2] and Ar = (2.6 ± 3.6 ± 
0.8) X 10-3 [8], respectively. They are consistent with zero, hence showing no indication 
of CP violation. 



2 Data selection 

LHCb is a precision heavy flavour experiment which exploits the abundance of charm par- 
ticles produced in collisions at the Large Hadron Collider (LHC). The LHCb detector [5] 
is a single arm spectrometer at the LHC with a pseudorapidity acceptance of 2 < < 5 
for charged particles. High precision measurements of flight distances are provided by the 
Vertex Locator (VELO), which consists of two halves with a series of semi-circular silicon 
microstrip detectors. The VELO measurements, together with momentum information 
from forward tracking stations and a 4 Tm dipole magnet, lead to decay-time resolutions 
of the order of one tenth of the lifetime. Two Ring-Imaging Cherenkov (RICH) de- 
tectors using three different radiators provide excellent pion-kaon separation over the full 
momentum range of interest. The detector is completed by hadronic and electromagnetic 
calorimeters and muon stations. The measurements presented here are based on a data 
sample corresponding to an integrated luminosity of 29 pb"^ of pp collisions at y/s = 7 TeV 
recorded during the LHC run in 2010. 

2.1 Trigger selection 

The LHCb trigger consists of hardware and software (HLT) stages. The hardware trigger 
is responsible for reducing the LHC pp interaction rate from (9(10) MHz to the rate at 
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which the LHCb subdetectors can be read out, nominally 1 MHz. It selects events based 
on the transverse momentum of track segments in the muon stations, the transverse energy 
of clusters in the calorimeters, and overall event multiplicity. 

The HLT further reduced the event rate to about 2 kHz in 2010, at which the data was 
stored for offline processing. The HLT runs the same software for the track reconstruction 
and event selection as is used offline and has access to the full event information. 

The first part of the HLT is based on the reconstruction of tracks and primary inter- 
action vertices in the VELO. Heavy flavour decays are identified by their large lifetimes, 
which cause their daughter tracks to be displaced from the primary interaction. The 
trigger first selects VELO tracks whose distance of closest approach to any primary in- 
teraction, known as the impact parameter (IP), exceeds 110 (xm. In addition the tracks 
are required to have at least ten hits in the VELO to reduce further the accepted rate 
of events. This cut limits the fiducial volume for D° decays and therefore rejects events 
where the candidate has a large transverse component of the distance of fiight, causing 
an upper bound on the decay-time acceptance. The term decay-time acceptance will be 
used throughout this paper to refer to the selection efficiency as a function of the D° 
decay time. Selected tracks are then used to define a region of interest in the tracking 
stations after the dipole magnet, whose size is defined by an assumed minimum track 
momentum of 8 GeV/ c; hits inside these search regions are used to form tracks traversing 
the full tracking system. Tracks passing this selection are fitted, yielding a full covariance 
matrix, and a final selection is made based on the track-fit quality and the track x^(IP). 
The x^(IP) is a measure of the consistency with the hypothesis that the IP is equal to 
zero. At least one good track is required for the event to be accepted. The requirements 
on both the track IP and on the x^(IP) reduce the number of candidates with a short 
decay time. 

In the second part of the HLT, an exclusive selection of candidates is performed 
by reconstructing two-track vertices. Further cuts are placed on the x^(IP) of the D° 
daughters and the displacement significance of the vertex from the primary interaction, 
as well as a requirement which limits the coUinearity angle between the momentum 
and the direction of fiight, as defined by the primary and decay vertices. These cuts 
all affect the distribution of the decay time of the candidates. Additional cuts are 
placed on track and vertex fit quality, and on kinematic quantities such as the transverse 
momentum of the candidate, which have no effect on the decay-time distribution. 

2.2 Offline selection 

Given the abundance of charm decays, the selection has been designed to achieve high 
purity. It uses similar requirements to those made in the trigger selection, though often 
with tighter thresholds. In addition it makes use of the RICH information for separating 
kaons and pions to achieve a low misidentification rate. A mass window of ±16MeV/c^ 
(about ±2(t) is applied to the invariant mass of the two daughter particles using the 
appropriate mass hypotheses. After these criteria have been applied there is negligible 
remaining cross-feed between the different two-body decay modes. 
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Flavour tagging of the decays is done by reconstructing the D*'^ — j- D^tt^ decay, 
where the charge of the slow pion, tTs, determines the flavour of the meson at produc- 
tion. The selection applies loose requirements on the kinematics of the bachelor pion and 
the quality of the D*^ vertex fit. The most powerful variable for selecting the D*"*" decay 
is the difference in the reconstructed invariant masses of the D*^ and the candidates, 
Am. Candidates are required to have Am in the range | Am — 145.4 MeV/c^| < 2.0MeV/c^. 

Events with multiple signal candidates are excluded from the analysis. For tagged 
decays this causes a reduction of the number of candidates of about 15% due to the high 
probability of assigning a random slow pion to form a D*~^ candidate. The numbers of 
selected candidates are 286,155 for )■ K'-k^ and 39,262 for K^K~ decays. 

3 Determination of proper-time acceptance effects 

Since absolute lifetime measurements are used to extract ycp and Ap, it is essential 
to correct for lifetime-biasing effects. The analysis uses a data-driven approach that 
calculates, for each candidate and at every possible decay time, an acceptance value 
of zero or one which is related to the trigger decision and offline selection. The final 
per-event acceptance function is used in the normalisation of the decay-time probability 
density function (PDF) as described in the following section. 

The method used to determine decay-time acceptance effects is based on the so-called 
"swimming" algorithm. This approach was first used at the NAll spectrometer [TUj . 
further developed within DELPHI [II] and CDF [I2l[l3], studied at LHCb [H US], and 
applied to the measurement of the -8° — )■ K^K~ lifetime [16] . 

The key to this method is the ability to execute the LHCb trigger software, including 
the reconstruction, in precisely the same configuration used during data taking. This is 
made possible by the implementation of all lifetime-biasing requirements of the trigger 
in software rather than in the hardware. The acceptance as a function of decay time is 
evaluated per event by artificially moving the position of the primary interaction vertices 
reconstructed in the trigger along the direction of the D*^ momentum in order to give the 
D*^ candidate a different decay time. In events containing multiple primary vertices, all 
are moved coherently in the direction of the D*^ momentum. An analogous procedure is 
used to obtain the decay-time acceptance of the offline selection. 

A decay-time acceptance function for any single event is in the simplest case a step 
function, as shown in Fig. [H since the kinematics and chosen decay time of the decay 
fully determine whether the event is triggered by this candidate or not. It is important 
to note that the acceptance function for a given event does not depend on the measured 
decay time of that event, tmeas- Accepted (rejected) regions take an acceptance value of 
1 (0). In this method decay-time independent selection efficiencies are factorised out and 
hence do not affect the result. The presence of additional interaction vertices can lead to 
regions of no acceptance and the VELO geometry puts an upper limit on the accepted 
range. Thus, a general decay-time acceptance function is given as a series of steps or top- 
hat functions. The decay times at which the event enters or leaves one of these top hats 
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are called turning points. The acceptance functions of the trigger and offline selections 
are combined to a single acceptance function by including only the ranges which have 
been accepted by both selections. 
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Figure 1: Variation of the decay-time acceptance function for a two-body decay when 
moving the primary vertex along the momentum vector. The shaded, light blue regions 
show the bands for accepting a track impact parameter. While the impact parameter of 
the negative track (IP2) is too low in (a) it reaches the accepted range in (b). The actual 
measured decay time, tmeas, lies in the accepted region which continues to larger decay 
times (c). 



The idea of studying the decay-time dependence of the acceptance in principle requires 
moving the hits produced by the decay products. The implementation of moving 
the primary vertices instead leads to significant technical simplifications. However, this 
procedure ignores the fact that events are no longer accepted if the mother particle has 
such a long decay time that one or both tracks can no longer be reconstructed inside the 
VELO. This is a very small effect as a meson has to fly ten to a hundred times its 
average distance of flight in order to escape detection in the VELO. Nevertheless, this 
effect can be estimated based on the knowledge of the position of the VELO modules and 
on the number of hits required to form a track. Using the information on the position 
of the VELO sensors, the limit of the acceptance is determined by swimming the tracks 
along the momentum vector. The result is treated as another per event decay-time 
acceptance and merged with the swimming results of the trigger and offline selections. 

Finally, the track reconstruction efficiency in the trigger is reduced compared to the 
offline reconstruction due to the requirements described in Sect. |2j It has been verified, 
using a smaller sample acquired without a lifetime biasing selection, that this relative 
reconstruction efflciency does not depend on the decay time of the candidate with a 
precision of 3 x 10~^, and therefore introduces no significant additional acceptance effect. 
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Figure 2: Am vs m^o distribution for — )■ K^n^ candidates. The contribution of ran- 
dom slow pions extends around the signal peak in the vertical direction while background 
is visible clS cL horizontal band. 

4 Fitting method 

The peak in Am from true D*~^ decays is parametrised as the sum of three Gaussians; 
two of which have a common mean and a third which has a slightly higher mean. The 
random tTs background PDF is given by 

, , f Am\^ ( , Am-d\ , /Am \ , . , . 
j^X^m)=[-—\ fl-exp( Am>rf, (6) 

where a and h define the slope at high values of Am, c defines the curvature at low 
values of Am and Am = d defines the threshold below which the function is equal to 
zero. Figure [2] shows the Am vs m^io distribution and Fig. E] shows the fit to the mass 
difference between the reconstructed invariant masses of D*'^ and D'^ candidates, Am. 

The signal yield is extracted from fits to the reconstructed invariant mass distri- 
bution after application of the cut in Am. The fit model for the signal peak has been 
chosen to be a double Gaussian and background is modelled as a first-order polynomial. 
The background level is evaluated to be about 1% for K'n'^ decays and about 

3% for — )■ K^K~ decays. It consists of combinatorial background and partially re- 
constructed or misidentified decays. If the latter stem from a D*^ decay they have a 
peaking distribution in Am similar to signal candidates. The data in the mass sidebands 
are insufficient to reliably describe the background shape in other variables, so the back- 
ground contribution is neglected in the time- dependent fit and a systematic uncertainty 
is estimated accordingly. 

Events inside the signal windows in Am and m^jo are used in the lifetime fit, where 
mesons produced at the primary vertex (prompt) have to be distinguished from those 
originating from h hadron decays (secondary). The combined PDF for this decay-time 



6 




140 145 150 155 160 140 145 150 155 160 

Am (MeV/c^) Am (MeV/c^) 



Figure 3: Am fit projections of (left) D'^ — )> K^ir^ and (right) D° — )■ K^K~ candidates. 
Shown are data (points), the total fit (green, solid) and the background component (blue, 
dot-dashed). 

dependent fit is factorized as 
/(x'(IPD),t,A)= /ip(x'(IPD)|t, A class) /^(tlA class) /Tp(A|class)P(class). (7) 

class 
=prompt, 
secondary 

The four factors on the right-hand side of Eq. [TJ which will be described in detail below, 
are: 

• the time-dependent PDFs for the lnx^(IP£)) values for prompt and secondary 
mesons; 

• the decay-time PDFs for prompt and secondary mesons; 

• the PDF for the turning points which define the acceptance A\ 

• the fractions of prompt and secondary decays among the signal candidates. 

The separation of prompt and secondary mesons is done on a statistical basis 
using the impact parameter of the candidate with respect to the primary vertex, IP^^. 
For prompt decays, this is zero up to resolution effects, but can acquire larger values 
for secondary decays as the candidate does not in general point back to the primary 
vertex. Given an estimate of the vertex resolution is available on an event-by-event basis, 
it is advantageous to use the of the IP o instead of the impact parameter value itself. 
The natural logarithm of this quantity, ln(x^(IP£))), allows for an easier parametrisation. 
Empirically, the sum of two bifurcated Gaussians, i.e. Gaussians with different widths 
on each side of the mean, and a third, symmetric Gaussian, all sharing a common peak 
position, is found to be a suitable model to describe the ln(x^(IP£))) distribution for both 
prompt and secondary . 
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For the prompt class the \n{x^ [IP £>)) distribution does not change with decay 
time as the true value is zero at all times and the resolution of IPd can be assumed to 
be independent of the measured decay time. For secondary decays the decay-time 
and ln(x^(IP£))) are correlated. The width of the ln(x^(IP£))) distribution is found to 
be approximately constant in decay time for both prompt and secondary mesons. 
As Monte Carlo simulation studies suggest that secondary decays have a larger width 
in this variable, a scale factor between the widths for prompt and secondary mesons is 
introduced. The mean value of ln(x^(IP/))) increases with decay time, which reflects 
the fact that mesons coming from other long-lived decays do not necessarily point 
back to the primary vertex and that they may point further away the further their parent 
particle flies. The functional form for this time dependence is based on simulation and 
all parameters are determined in the fit to data. 

The decay-time PDF, ft{t\A, class) is modelled as a single exponential for the prompt 

class and as a convolution of two exponentials for secondary decays. To account for 
resolution effects, these are convolved with a single Gaussian resolution function. The 
parameters of the resolution model are obtained from a fit to the decay time distribution 
of prompt J/ip events. The resulting dilution is equivalent to that of a single Gaussian 
with a width of 50 fs [17] . The decay-time probability densities are properly normalized 
by integrating their product with the acceptance function A, evaluated by the swimming 
method, only over the decay-time intervals for which the event would have been accepted. 
Hence, the acceptance turning points are used as boundaries in the integration. 

Finally, a PDF for the per-event acceptance function is needed. While the first ac- 
ceptance turning point, i.e. the one with the smallest decay time, depends on the 
decay topology, the others are governed more by the underlying event structure, e.g. 
the distribution of primary vertices. The primary vertex distribution is independent of 
whether the candidate is of prompt or secondary origin. Hence, the PDF can be ap- 
proximated as /Tp(A|class) ~ /tp (TP i| class), where TPi denotes the position of the first 
turning point. The distribution for /Tp(TPi|prompt) is obtained by applying a cut at 
In (IPd) < 1, thus selecting a very pure sample of prompt decays. The distribution for 
/Tp(TPi|secondary) is obtained from the distribution of TPi weighted by the probability 
of each candidate being of secondary decay origin. 

An initial fit is performed using the full In {IPd) distribution and all parameters in 
the description of this term are then fixed in the final fit. A cut is then applied requiring 
In x^ (IPd) < 2 in order to suppress the fractions of both background and secondary 
candidates to less than a few percent. The final fit is performed on this reduced sample. 
The effect of this procedure is estimated in the systematic uncertainty evaluation. 

5 Cross-checks and systematic uncertainties 

The method for absolute lifetime measurements described in Sect. H] comprises three main 
parts whose accuracy and potential for biasing the measurement have to be evaluated in 
detail: 
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• the determination of the event-by-event decay-time acceptance; 

• the separation of prompt from secondary charm decays; 

• the estimation of the decay time distribution of combinatorial background. 

Since the contribution of combinatorial background is ignored in the fit, it is important to 
evaluate the corresponding systematic uncertainty. Furthermore, several other parameters 
are used in the fit whose systematic effects have to be evaluated, e.g. the description of 
the decay-time resolution. It is generally expected that the systematic uncertainties in 
ycp are similar to or larger than those in as in ycp two different final states contribute 
to the measurement. 

Several consistency checks are performed by splitting the dataset into subsets. The 
stability is tested as a function of run period, D° momentum and transverse momen- 
tum, and primary vertex multiplicity. No significant trend is observed and therefore no 
systematic uncertainty assigned. 

The fitting procedure is verified using simplified Monte Carlo simulation studies. No 
indication of a bias is observed and the statistical uncertainties are estimated accurately. 

As an additional check, a control measurement is performed using the lifetime asym- 
metry of and decays to the Cabibbo favoured decay )■ K~7r~^. The result is in 
agreement with zero and the flavour-averaged lifetime is found to be consistent with 
the world average. Detailed results are given in Sect. [Hi The fit results for )■ K^K~ 
decays were not looked at throughout the development of the method and the study of 
systematic uncertainties for the analyses of ycp and A^. 

5.1 Evaluation of systematic uncertainties 

Particle decay times are measured from the distance between the primary vertex and 
secondary decay vertex in the VELO. The systematic uncertainty from the distance scale 
is determined by considering the potential error on the length scale of the detector from 
the mechanical survey, thermal expansion and the current alignment precision. A relative 
systematic uncertainty of 0.1% is assigned to the measurements of absolute lifetimes, 
translating into a relative uncertainty of 0.1% on Ay and ycp- 

The method to evaluate the turning points of the decay-time acceptance functions 
described in Sect. |3] uses an iterative approach which estimates the turning points to a 
precision of about 1 fs. Two scenarios have been tested: a common bias of all acceptance 
turning points and a common length scaling of the turning points, which could originate 
from differences in the length scale in the trigger and offline reconstructions. From a 
variation of the bias and the scale, a systematic uncertainty of 0.1 x 10~^ on Ay and ycp 
is determined. 

The reconstruction acceptance is dominated by the VELO geometry, which is ac- 
counted for by the method described in Sect. [31 This leads to a correction of less than 1 fs 
on the absolute lifetime measurements, i.e. a relative correction of about 0.24%. No fur- 
ther systematic uncertainty is assigned to Ay or ycp as the size of this relative correction 
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is negligible. Additional studies of the reconstruction efficiency as a function of variables 
governing the decay geometry did not provide any indication of lifetime biasing effects. 

The decay-time resolution is modelled by a single Gaussian. The width of the resolu- 
tion function is varied from its nominal value of 0.05 ps between 0.03 ps and 0.07 ps. The 
range of variation was chosen to cover possible alignment effects as well as effects from 
the different final state used to evaluate the resolution. The result leads to a systematic 
uncertainty of 0.1 x 10~^ for Ar and ycp- 

The fit range in decay time is restricted by lower and upper limits. The lower limit is 
put in place to avoid instabilities in regions with extremely low decay-time acceptances 
and very few events. The default cut value is 0.25 ps which is close to the lower end of the 
observed range of events. This cut is varied to both 0.2 ps and 0.3 ps. The result leads to 
a systematic uncertainty of 0.1 x 10"'^ for A-p and 0.8 x 10~^ for ycp- 

The upper limit of the ffi range in decay time is put in place to minimise the impact of 
long-lived background events. The default cut is put at 6 ps which corresponds to about 
15 lifetimes. This cut is varied to 5ps and 8ps. The result leads to a systematic 
uncertainty of 0.2 x 10~^ for A-p and ycp- 

The description of the contribution from combinatorial background is studied by vary- 
ing its relative amount in the data sample and repeating the ffi. This is done by changing 
the Am window from the default of ±2 MeV/c^ to ±1 MeV/c^ and ±3 MeV/c^. The result 
leads to a systematic uncertainty of 1.3 x 10~^ for A^ and 0.8 x 10"'^ for ycp- 

Events that originate from secondary charm decays are the background with the largest 
impact on the fit procedure as they have a very different decay-time distribution compared 
to prompt charm decays, but they peak in the invariant mass and Am distributions. 
Also a fraction of combinatorial background events appear to be secondary-like in their 
lnx^(IP£)) distribution. The cut of lnx^(IP/)) < 2 removes a large fraction of secondary- 
like events. However, it is important that the remainder is properly modelled and does 
not bias the signal lifetime. Varying this cut changes the relative number of secondary- like 
decays in the sample and therefore tests the stability of the secondary description in the 
fit model. The fraction of secondary-like combinatorial background events is also altered 
with this test. The lnx^(IP£)) cut is varied from 1.5 which is just above the peak of 
the prompt distribution to 3.5 where the probability densities for prompt and secondary 
decays are about equal. The result leads to a systematic uncertainty of 1.6 x 10~^ for A^ 
and 3.9 x 10~^ for ycp- The uncertainty is significantly larger for ycp than for Ar as may 
be expected from the difference in the background level in the channels involved in the 
ycp measurement. 

Additional studies were performed to estimate the potential impact of neglecting back- 
ground events in the ffi. A background component was added to a simplified simulation. 
The background decay time distribution was generated using extreme values of fits to the 
distribution observed in mass sidebands. The average bias on the measurement of ycp 
was about 2 x 10~^. Since this is consistent with the assigned systematic uncertainty, we 
do not assign any additional uncertainty. 

Furthermore, a background component was added to the decay-time PDF with a 
fixed fraction and average lifetime. The fraction of this component, which was assumed 
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Table 1: Summary of systematic uncertainties. 



Effect 


At (10-3) 


ycp (10-=^) 


Decay-time acceptance correction 


0.1 


0.1 


Decay-time resolution 


0.1 


0.1 


Minimum decay-time cut 


0.1 


0.8 


Maximum decay-time cut 


0.2 


0.2 


Combinatorial background 


1.3 


0.8 


Secondary-like background 


1.6 


3.9 


Total 


2.1 


4.1 



to be secondary-like, was varied. A change in the fit result for ycp of (all background 
secondary-like) to 4 x 10"^ (all background prompt-like) was observed. As it is known that 
a fraction of the background events are secondary-like, this result is considered consistent 
with the simplified simulation results. 

5.2 Summary of systematic uncertainties 

Table [1] summarises the systematic uncertainties evaluated as described above. The main 
systematic uncertainties are due to neglecting the combinatorial background and to the 
contribution of secondary-like decays. The total systematic uncertainties for At and ycp-, 
obtained by combining all sources in quadrature, are 2.1 x 10"^ and 4.1 x 10"^, respectively. 

6 Results and conclusion 

The measurement of ycp is based on absolute lifetime measurements as described in 
Sect, m It uses flavour-tagged events reconstructed in the decay chain D*^ D^n^, 
with D° and D° decays fitted simultaneously per decay mode. The Inx^ilPn) projection 
of the final fit is shown in Fig. HI 

The result for the lifetime measured in ^ K^n^ decays is t{D^) = 410.2 ± 0.9 fs 
where the uncertainty is statistical only. The result for the lifetime is found to be in 
agreement with the current world average [IB]. Combining with the ^ K^K~ lifetime 
measurement, t{D^) = 408.0 ± 2.4statfs, this leads to the final result for ycp of 

l/cp= (5.5 ±6.3stat± 4. Isyst) X 10-\ 

The measurement of At is performed based on the same dataset and applying the same 
fitting method as used for the extraction of ycp- A control measurement is performed 
using decays to the Cabibbo favoured mode Z)"— >■ K~7r~^ by forming a lifetime asymmetry 
analogous to Eq. |2j The measured flavour-tagged lifetimes are effective parameters since 
the fitted distributions also include mistagged events. For the control measurement using 

—7- K~7T~^ decays this contamination is ignored as it is negligible due to the Cabibbo 
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Figure 4: lnx^{lPD) fit projection of — )■ K^K~ candidates in logarithmic scale. Shown 
are data (points), the total fit (green, solid), the prompt signal (blue, short-dashed), and 
the secondary signal (purple, long-dashed). 

suppression of the mistagged decays. The result for the asymmetry is A^^''^^ = (—0.9 ± 
2.2stat) X 10~^ which is consistent with zero, according to expectations. 

For the extraction of Ar, the mistagged decays are taken into account by expressing 
the measured effective lifetimes, t'^^, in terms of the flavour-tagged lifetimes, r(il'°) and 
t{D^), and the mistag rate, e±, where the sign is according to the sign of the tagging 
pion: 



The mistag rates are assumed to be independent of the flnal state and are extracted from 
the favoured )■ K~'ir~^ decays as half the fraction of the random slow pion background 
in the signal region of the Am distribution. They are found to be about 1.8%. The 
systematic uncertainty due to this correction is negligible. 

The projection of the decay-time flt to and D'^ candidates in — )■ K^K~ decays 
is shown in Fig. O After applying the mistag correction, the resulting value of Ay is 



Both results on ycp and Ay are compatible with zero and in agreement with previous 
measurements [21 13 E]- Future updates are expected to lead to signiflcant improvements 
in the sensitivity. The systematic uncertainty is expected to be reduced by an improved 
treatment of background events which will be possible for the data taken in 2011. 
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